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This paper develops a simplified set of models describing asexual and sexual replication in unicel- 
lular diploid organisms. The models assume organisms whose genomes consist of two chromosomes, 
where each chromosome is assumed to be functional if it is equal to some master sequence oq, and 
non-functional otherwise. The first-order growth rate constant, or fitness, of an organism, is de- 
termined by whether it has zero, one, or two functional chromosomes in its genome. We assume 
that an organism with no functional chromosomes has zero fitness. For a population replicating 
asexually, a given cell replicates both of its chromosomes, and then the cell divides and splits its 
genetic material evenly between the two cells. For a population replicating sexually, a given cell 
first divides into two haploids, which enter a haploid pool. Within the haploid pool, haploids fuse 
into diploids, which then divide via the normal mitotic process. Haploid fusion is modeled as a 
second-order rate process. We review the previously studied case of selective mating, where it is 
assumed that only haploids with functional chromosomes can fuse, and also consider the case of ran- 
dom haploid fusion. When the cost for sex is small, as measured by the ratio of the characteristic 
haploid fusion time to the characteristic growth time, we find that sexual replication with random 
haploid fusion leads to a greater mean fitness for the population than a purely asexual strategy. 
However, independently of the cost for sex, we find that sexual replication with a selective mating 
strategy leads to a higher mean fitness than the random mating strategy. This result is based on 
the assumption that a selective mating strategy does not have any additional time or energy costs 
over the random mating strategy, an assumption that is discussed in the paper. The results of this 
paper are consistent with previous studies suggesting that sex is favored at intermediate mutation 
rates, for slowly replicating organisms, and at high population densities. 
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I. INTRODUCTION 

The evolution and maintenance of sexual replication 
is one of the central questions in modern evolutionary 
biology [TJ [21 El IU [5] . The theories with the broadest ac- 
ceptance are that sex allows for the removal of deleterious 
mutations from a population [HlHllI], or that sex allows 
for faster adaptation in dynamic environments [1] [8] . 

As has been noted in previous studies [HI [111], the var- 
ious theories for the existence of sex are incomplete, in 
that they do not explain why some organisms are obli- 
gatory sexual, while other organisms cither alternate be- 
tween asexual and sexual replication, or are asexual repli- 
cators with some ability for recombination with other or- 
ganisms. 

In a recent set of papers [HI [101 HU H], Tannen- 
baum, Fontanari, and Lee studied the competition be- 
tween asexual and sexual replicators in various regimes. 
Two of the models consider single-celled organisms, that 
replicate by dividing into haploids. These haploids then 
enter a haploid pool, where they fuse with other haploids. 
The resulting diploid then divides through the normal 
mitotic process. 

Two other models were developed to be more appropri- 
ate for modeling multicellular organisms, which release 
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either asexual spores or gametes. Here, we considered the 
case where the population produces identical gametes, 
and where the population produces distinct sperm and 
egg gametes. 

For both sets of models it was found that sexual repli- 
cation is favored in slowly replicating organisms and high 
population densities. For the case of multicellular organ- 
isms, it was found that distinct sperm and egg gametes 
were necessary to maintain a selective advantage for sex- 
ual replication over asexual replication |12j . 

In all the models considered, the authors assumed that 
there is a time cost for sex, due to the time it takes for 
a given haploid to find another haploid with which to 
fuse. When the time cost for sex is low, then the selective 
advantage for sex outweighs the fitness penalty. Once the 
time cost for sex becomes sufficiently large, the selective 
advantage for sex no longer outweighs the fitness penalty, 
and asexual replication becomes the preferred replication 
strategy. 

In order to facilitate an exact analysis of the models, 
the authors made a number of simplifiying assumptions: 
It was assumed that the organisms have diploid genomes 
consisting of only two chromosomes, and that each chro- 
mosome is only considered functional, or viable, if it is 
equal to some master sequence. It was also assumed that 
the fitness of the organisms is determined by the number 
of viable chromosomes in the genome (0, 1, 2), and that 
organisms with two defective chromosomes have zero fit- 
ness. 
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Finally, one of the key simplifications that was made 
in the sexual replication models was the assumption that 
only haploids with viable chromosomes can fuse with 
one another. The reasoning behind this selective mating 
strategy was that haploids with defective chromosomes 
are not viable, and simply cannot participate further in 
the replication process. 

While the assumption of a selective mating strategy is 
the easiest one to study, and while it is broadly consistent 
with the observation that organisms do engage in mate 
selection, it is nevertheless an overly restrictive assump- 
tion that must be relaxed if one wants to develop more 
realistic sexual replication models. The reason for this is 
that it is impossible, in practice, for one organism to read 
another organism's genome. Even if this were possible, 
the correlation between the genome sequence and fitness 
is extremely difficult to obtain. 

Organisms that engage in mate selection look for cer- 
tain behaviors and physical attributes, known as indi- 
cator traits, that suggest that the given organism has a 
good genome |13j . However, the map from such traits 
to genome is not exact. Furthermore, the more accuracy 
one wants in assessing the fitness of an organism, the 
more time one has to spend studying the organism. For, 
if assessing the fitness of an organism is equivalent to 
reading its genome, then the more accuracy one wants, 
the more of the genome that has to be read. 

Therefore, a selective mating strategy has an addi- 
tional time cost over other mating strategies, leading to a 
fitness penalty that may eliminate the advantage for sex- 
ual replication entirely. This time cost was not explicitly 
considered in previous studies. 

Even if the time cost for a selective mating strategy 
were explicitly considered, we have shown that there will 
always be an uncertainly as to the exact fitness of the 
organism. As a result, more realistic models for sexual 
replication need to consider mating strategies that take 
into account the uncertainty that one organism has of an- 
other's fitness. Otherwise, the conclusions that we have 
drawn regarding the selective advantage for sex are open 
to the criticism that they are based on an unrealistic and 
overly restrictive mating strategy. 

In this paper, we take the opposite extreme from a se- 
lective mating strategy, and consider sexual replication 
with a random mating strategy. That is, we assume that 
all haploids participate in the replication process, and 
fuse with one another at random. As long as the cost for 
sex is negligible, we find that, even with this non-selective 
mating strategy, the mean fitness for a sexually replicat- 
ing population is greater than that of the correspond- 
ing asexually replicating population. Nevertheless, if the 
time cost associated with the selective mating strategy 
is negligible as well, then the mean fitness of the selec- 
tive mating strategy is greater than that of the random 
mating strategy, independent of the cost for sex. As the 
cost for sex increases, sexual replication via either the 
selective or random mating strategies only outcompetes 
asexual replication over progressively smaller ranges of 



replication fidelities, and ceases to be advantageous en- 
tirely once the cost for sex crosses a threshold value. 

This paper is organized as follows: In the following sec- 
tion (Section II), we develop and analyze the mutation- 
selection equations appropriate for analyzing an asexu- 
ally replicating, two-chromosomcd, diploid population. 
In Section III, we develop the mutation-selection equa- 
tions appropriate for analyzing the sexually replicating 
analogue of the population considered in Section II. We 
consider both selective and random mating strategies. 
Although the selective mating strategy has already been 
studied, we review it here for the sake of completeness 
and consistency of notation. In Section IV, we compare 
all three replication strategies, and show that a random 
mating strategy outcompetes asexual replication when 
the cost for sex is negligible. We also show that when 
the cost of implementing the selective mating strategy 
is negligible, then the selective mating strategy outcom- 
petes the random mating strategy, independently of the 
cost for sex. Section V is the Discussion, where we re- 
view and discuss the implications of our results. Finally, 
in Section VI we conclude the paper with a summary of 
the main results of the paper, as well as our plans for 
future research. 



II. THE ASEXUAL REPLICATION MODEL 

We consider a population of unicellular organisms 
whose genomes consist of two chromosomes. We assume 
that a given chromosome, denoted <r, is functional if and 
only if it is equal to some "master" , or wild- type, chro- 
mosome (To- Assuming first-order exponential growth, we 
then assume that the fitness, or first-order growth rate 
constant, of a given genome is determined by whether 
the genome consists of zero, one, or two functional chro- 
mosomes. To this end, we let k vv denote the first-order 
growth rate constant of organisms with two functional, 
or equivalently, viable chromosomes. We let k vu denote 
the first-order growth rate constant of organisms with 
one functional and one non-functional (unviable) chro- 
mosome. Finally, we let n uu denote the first-order growth 
rate constant of organisms with two non-functional chro- 
mosomes. 

We assume that n vv > n vu > n uu . We will also assume 
that k uu — 0, which makes sense, since an organism with 
two defective chromosomes is not expected to grow. We 
also define k = K vv , and a — k vu /k vv . 

We now divide the populations into three subpopu- 
lations: We let n VVl n VUl and n uu denote the number 
of organisms with two, one, and zero viable chromo- 
somes. We let n be the total number of organisms, so 
that n = n vv + n vu + n uu . The population fractions 
x vv , x vu , and x uu are then defined via x vv — n vv /n, 
Xvu — ^i vu /n : and x uu — ti uu I ji. 

To develop the mutation-selection equations describing 
the evolutionary dynamics of the asexually replicating 
equations, we assume that the replication of each cell 
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occurs as follows: The two chromosomes line up along the 
center of the cell and replicate (see Figure 1). Because 
replication is in general not error-free, we let p denote the 
probability that a given genome is replicated correctly. If 
the genome is sufficiently long, then the probability that 
a mutation will occur in a previously mutated region of 
the genome is negligible, so that an unviable chromosome 
produces an unviable daughter with probability 1. This 
assumption is known as the neglect of backmutations. 

We also let s denote a co-segregation parameter, which 
is simply the probability that a parent chromosome co- 
segregates with the other parent chromosome in the cell. 
Figure 1 illustrates the various parent cell configurations 
and the final daughter cell configurations, along with 
their associated probabilities. It should be noted that 
1 — s is the same as the r parameter defined in [10] . 

With these definitions in hand, we may develop ex- 
pressions for dn vv /dt 7 dn vu /dt, and dn uu /dt. Changing 
variables from population numbers to population frac- 
tions, we obtain the mutation-selection equations, 



dt 

dt 

dx uu 



dt 1 

+ S(1 - p) 2 K vv X, 



(1) 



where /^(^) = ('i-/Ti)(dTi/dtj — Kw-^w ~\~ Kvu^vu ^uu^v 
and, 



A vv = l -[2p + s{l-p) 2 ] 
B vv = [l-p][l-s{l-p)] 



(2) 



The quantity R(t) is the mean fitness of the population, 
since it measures the first-order growth rate of the popu- 
lation as a whole. In order to determine which replication 
strategy is advantageous in a given regime, we compare 
the steady-state mean fitnesses of the populations em- 
ploying the various strategies. The population with the 
largest mean fitness will drive the others to extinction, 
and so the corresponding replication strategy is the ad- 
vantageous one. 

From quasispecies theory, it may be shown that the 
above system of equations converges to a steady-state, 
and that the steady-state mean fitness is given by R(t = 
oo) = max{(2A vv — 1)k vv , K vu p}, assuming that k uu — 
[14] . This implies that there exists a p cr it € [0, 1] such 
that R(t = oo) = 2A VV — 1 for p £ [p C rit, 1], and R{t = 

OO) = K vu p for p e [0, Pcrit)- 

If we define <p a = K,(t = oo)/k vv , then we have that 
4> a = 2A VV - 1 for p e [pcrit, 1], and ap for p g [0,p C rit], 
where p cr it is defined by the equality ap — 2A VV — 1. 

Note that (f> a is simply the steady-state mean fitness 
of the asexual population, normalized by the fitness of 



the wild-type, k vv . When we analyze the sexual replica- 
tion models, we will also work with the normalized mean 
fitnesses, as it will prove convenient to do so. When 
comparing the various replication strategies, the strategy 
with the largest normalized mean fitness at steady-state 
is the one that will outcompete the others for the given 
set of parameters. 



III. THE SEXUAL REPLICATION MODELS 
A. The general model 

The general sexual replication model we are consider- 
ing is illustrated in Figure 2: A diploid cell grows to ma- 
ture size, with a first-order growth rate constant given by 
Kw, or k uu — 0, depending on whether the genome 
has two, one, or zero functional chromosomes, respec- 
tively. The mature diploid then divides into two hap- 
loids, which enter a haploid pool. The haploids fuse with 
one another, and the resulting diploids then immediately 
divide via the normal mitotic process. 

It should be noted that the corresponding first-order 
growth rate constants for the asexual and sexual popula- 
tions are taken to be equal, since the first-order growth 
rate constants measure the characteristic time it takes a 
diploid to double in size before dividing. 

With sexual replication, it is necessary to keep track of 
the haploid as well as the diploid population. Therefore, 
in addition to the quantities n vv , n vu , and n uu defined in 
the previous section, we also have the quantities n v and 
n u , corresponding to the number of viable and unviable 
haploids, respectively. 

The haploid fusion process is modeled as a binary col- 
lision reaction characterized by second-order rate con- 
stants jvv, 7««j and j uu , corresponding to the v — v, v — u, 
and u — u haploid collisions. It should be noted that j vv 
and 7„ u are defined so that ( / y vv /V)ny and (7uu/V)n£ 
are the rates of disappearance of the v haploids and u 
haploids respectively, due to v — v and u — u haploid 
fusion respectively. The quantity j V u is defined so that 
(lvu/V) n vn u is the rate of disappearance of either the v 
haploids or the u haploids due to v — u haploid fusion. 

We also assume that the system volume increases so as 
to maintain a constant density of genomes in the popula- 
tion. That is, if V denotes the volume of the system, then 
we assume that p = [n vv + n vu + n uu + (1 /2)(n v + n u )]/V 
is constant. 

We now make the following definitions: We define 
n = n vv + n vu + n uu to be the total population of 
diploids, and R(t) = (l/n)(dn/dt) to be the mean fit- 
ness of the diploid population. We define the popula- 
tion ratios X V v — ^vv/^i ^vu — ^vu/^i ^uu — tl'uu/^i 

x v = n v /n and x u = n u /n. Finally, we define a diploid 
density p* = n/V = p/(l + (l/2)(x v + x u ). 

If we write down the differential equations governing 
the values of n vv , n vu , n uu , n v and n u , then changing 
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variables to the population ratios gives the mutation- 
selection equations, 



d>X V y 

dt 

dXyu 

dt 

dx uu , ~(4\\ _L * C 1 2 

— — T" )X UU ~\- 'JwP ^VV^y 

lit 

-% 

dx u 
~dt 



f\{t)X V -\~ QlKyyXyy ~\~ KyyXy^ 

IfvvP x v ~ ~~)vuP x v^u 

= -R{t)x u + 



^JvuP %v%u ^JuuP 



(3) 



where 



KyyXyy ^VU^VU ^UU^UU 

«) (4) 

and where we have defined the additional quantities , 
C vv , and Cyu via, 



1 



P 



1 



P 



(5) 



Although we are defining the mean fitness of the pop- 
ulation with respect to the diploid organisms, at steady- 
state any two mean fitnesses defined with respect to two 
distinct sets of subpopulations will be equal. The reason 
for this is that, at steady-state, the population reaches a 
mutation-selection balance, so that any two given sets of 
subpopulations will be in a fixed proportion to one an- 
other (e.g. the total diploid population versus the total 
population). Therefore, the per capita rate of increase 
of one subpopulation is equal to the per capita rate of 
increase of another. 

As a result, there is no ambiguity as to which sexual 
mean fitness to use when comparing whether a given sex- 
ual population outcompetes another sexual population, 
or an asexual population. 



B. Steady-state for the selective mating strategy 

The selective mating strategy is defined by 7„„ = 7, 
Ivu = luu = 0. This implies that only the v haploids 
are allowed to mate, while the u haploids are essentially 
thrown away. One justification for this mating strategy 
is that the u haploids contain defective genomes, and so 
are simply physically incapable of participating further in 
the replication process. Another justification, one that is 



likely more relevant to actual organisms, is that the v 
haploids have a way of determining the fitness of a po- 
tential haploid mate (via "indicator" traits, for example), 
and choose to only mate with the fittest haploids. 

For the selective mating strategy, the steady-state 
equations are given by, 

„*™2 



= -{k vv + n(t = oo))x vv + jp*x v A vv 
= -(«;„« + K(t = oo))x vu + -)p*x 2 v B vv 

K,(t OO^jXy ~\- *2k, vv X vv ~\- fcyyXyn 

K,(t OO^jX u ~\- fcyu^VU 



* 2 

IP x v 



(6) 



where n{t = 00) 



rv vv X v 



+ jp*x„. We pur- 



posely neglect the steady-state equation corresponding 
to dx uu /dt — 0, since this equation will not be necessary 
to determine the mean fitness of the population. 

Solving the last equation for jp*x% and substituting 
the result into the expression for R(t = 00) gives, 



OO^ KiyyXyy t%{t OC)) X y 



and so, 



x v 



K VV X V 

R(t = 00) 



- 1 



We also have, 



R(t = 00) 



(7) 



(8) 



(9) 



Substituting these values into the expression for R(t 
00), and making use of the fact that p* = p/(l 
(1/2) (x v + x u )) gives, after some manipulation, 



(x vv <fiss) 



<j) ss {x vv + ax vu + 4>Y 



1 K vv 

2 TP 



(10) 



where (j) ss is the normalized steady-state mean fitness of 
the population, and is given by (j) ss = R{t — oo)/k vv . 
Here, ss stands for "selective sexual". 

The steady-state equations for x vv and x vu then give, 



(1 + <p ss )x 
(a + (f) ss )x 



vv - 
vu 



■Ayy {^X 



vv ~r CW^vu 

yy ~\- OtXy U ~\~ 



+ 4>ss) 



Therefore, x vu /x v 
so that 

(1 ~t~ 4 > ss)x vv = A.', 



) (11) 

B vv /A vv X (1 + <f> aa )/(a + <j) sa ), 

, + (A vv + (1 + <t> 8S )B v 



a + < 



_ A vv (j) ss (a + cj) ss ) 

(1 + <f> as )(a + 4> ss - B vv a) - A vv (a + <f) ss ) 

Now, the first steady-state equation for x vv may be 
solved to give, x vv + ax vu + 4> ss = (1 + (j) ss )x vv / A vv , 
which may be substituted into Eq. (10) to give, 



" *)Xvv 

(12) 



A vv (l 



1 K r , 



5 (1 + SS ) 2 2 1P 



(13) 
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Substituting in the value for x vv gives, after some ma- 
nipulation, that, 



The steady-state equations for x vv and x vu give, 



(2A VV (1 - a) - 1 + ap)4>ss - ap] 2 _ 1 k vv 
<p ss (cl> ss + l) 2 (<p ss + a) 2 ~ 2^p~ 



(14) 



(1 $rs)Xvv — A vv (^ 



2x V v ~\~ OtXvu 

2(^X VV ~\- (XXyy^ 



) (0rs ~t~ X V v "t" Q^vu) 



We should note that the results we have obtained here 
for the selective mating strategy were previously derived 
by Tannenbaum and Fontanari [101 . Nevertheless, we 
reviewed them here for the sake of completeness. 



C. Steady-state for the random mating strategy 

The random mating strategy is defined by j vv = j vu = 
luu = 7- Here, any haploid pair is equally likely to fuse 
as any another haploid pair. 

The steady-state equations are, 

= -(Km, + R(t = oo))x vv + jp*x^A vv 

= -(k vu + R(t = oo))x vu + jp*xlB vv + -fp*x v x u B vu 

K.{t OO^j^Xy ~\- Xllj ~T~ ^(^KyyXyy ~\~ ^VU^VIl) 

--ip*(x v + x u ) 2 

K,(t OO^jXy -j - ^iKyyXyy ~\~ ^yyXvU TP X y{^X y ~f~ X y^j 

(15) 

where R(t = oo) = -k vv x vv - k vu x vu + r yp*(x v + x u ) 2 . 
The third equation is obtained by adding the equations 
corresponding to dx v /dt = dx u /dt = 0. 
We then have, 

R(t = oo) = k vv x vv + k vu x vu — R(t — oo)(x v + x u ) (16) 

so that, 

K VV X VV -\- fcvuXvu 



yCt -)- (prs )Xvu — J^vv \ / \ ) \(Prs < X V v < (%X VU ) 



+B V 



'Z^Xyy H - OtXyiij 
(^2iXyy ~\~ CXXyiijCXXy U 



A.{Xyy ~\~ Q.X VU ) 



Now, from Eq. (16) it can be seen that (j) rs < x v 
ax vu , and so from Eq. (20) we have, 



•&vv ~r OiXyu (prs I ^vv ^ 

V 2 7 P 

=r* Xyy H - OiXyu 0rS /(<Prs, A) 



(22) 



w here A = K vv /(2jp), /(0 rs ,A) = (1 + \f\(j> rs )/(l 
We now have, 



(1 + <t> rs )X vv = + /(0 rs , A)) 



(gOT, + 4>rsf(4>rs, A))" 

4</) rs /(</) rs ,A) 2 



/ i J, \ "Prs/C^rs, A) %vv j-, /-, ,,/ , .n 
(a + rs ) = -B™(1 + /(<prs, A)) X 

(x„„ + (jjrsfi^rs, A)) 2 



^r 



«(i = oo) 



1 



(17) 



4(pW(<pV s ,A) 2 

+B««(l + /(0r.,A)) x 
(</>rs/((?Vs, A) + Zm-)(0rs/((/>rs, A) - 
^rsf(4>rs, A) 2 

The second equation gives, 

(a + 0r S )(prs/((/>r S , A) - (a + cj) rs )x vv = 



(23) 



a 



l + /(0r„A) 



If we define x v — x v /(x v + x u ), x u = x u /(x v + x u ), 
then, 

and so, 



^^KyyXyy ~\~ f^VU^Vu) 
^(^tiyyXyy ~\~ ^VU^VIl) 



(19) 



By substituting the expression for X y | 3/ into the 
expression defining R(t = oo), we obtain, after some ma- 
nipulation, 



(x vv ~\~ OLX vu 0rs) 1 K vv 

4>rs(Xvv + ax vu + 4> rs ) 2 2 -fp 



(20) 



where 4> rs = n(t = oo)/n vv . Here, rs stands for "random 
sexual" . 



40rs/(0rs, A) 2 

[B vv ((f> rs f{(t) rs , A) + a;^,,) 2 + B vu ((j) 2 s f((j) rs ,\) 2 - x 2 vv )] 

(24) 

Now, 4> 2 s f((t> rSl \) 2 - Z 2 ,, = -(<p rs f(<f> rs ,X) + X vv ) 2 + 

2<t> 2 s f(4> rsi A) 2 + 20rs/(0r6, A)^^^, and so, after some ma- 
nipulation, we obtain, 

Xvv = 4>rsf(4>rs, A) X 

[(a + <j) rs )f(<j) rsi A) - l -aB V u{l + f(<p rs , A))]/ 

A)(l- 2a) 
+a(^S TO (l + /(0rs, A)) - /(^.., A))] 

(25) 

where we made use of the fact that B vv — 
B vu = -2A VV and 2(1 + (j> rs )xv V = (1 + 

f(<j>rs, \))/(2<j) rs f(<t> rs ,\) 2 )Avv(<f>rsf((i)rs,>*) + X V v) 2 ■ 
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Plugging the value of x vv back into the first equation 
from Eq. (23) we obtain, after tedious algebra, 

[l + <t> ra ] x 

[0rs/Or S , A) - a(-B vu (l + f((f> rs , A)) - f((f) rs ,X))] X 

A)(l - 2a) + a{^B vu {\ + f{cj> rs ,\)) - f(4> r , 
= A vv (l - a) 2 f(<f> rs , A)(l + f(<f>„, \))4> 2 rs 

Now, l+f((j) rs , A) = 2/ (1— \f\<f> rs ), and so, multiplying 
both sides by (1 — ^\cf) rs ) 2 gives, 

KiMrs + k 2 \f>47s- A 3 = (27) 

where, 

Ai = (<j) rs + l){<j) rs + a)((f) rs (l - 2a) - a) 

A 2 = 2(0 rs + 1)((1 - 2a)4> 2 rs - a 2 (l - p)4> rs + a 2 p) 

-2A VV {1 -a) 2 4> 2 rs 

A 3 = 2A vv (l- a ) 2 ct> 2 rs 

-{<t> rs + l)(0 rs - ap)(^rs(l - 2a) + ap) 

(28) 

The solution that is chosen is the one that gives A — ► 
as A3 — > 0, since A3 = is the equation defining the 
steady-state mean fitness when there is no cost for sex. 



IV. COMPARISON OF THE ASEXUAL AND 
SEXUAL REPLICATION STRATEGIES 

We now compare the various replication strategies. We 
consider first the case where there is no cost for sex, so 
that n vv /(^p) = 0, followed by the case where there is a 
non-zero cost for sex, so that n vv /("fp) > 0. 

The ratio K vv /(-fp) measures the cost for sex because 
it may be interpreted as the ratio of the characteristic 
time a haploid spends looking for another haploid with 
which to fuse, which is on the order of 1/(7/9), to the 
characteristic time it takes newly formed diploid cell to 
grow to maturity and divide, which is on the order of 
1/k vv . When this ratio is small, then the fraction of the 
organism's life cycle that is devoted to the haploid fusion 
process is small, so that the time cost associated with 
sex is small. Conversely, when this ratio is large, then 
the time cost associated with sex is large as well. 



A. Case 1: n vv /(-yp) = 

When n vv /("fp) = 0, the normalized mean fitness for 
the population replicating with the selective sexual repli- 
cation strategy is given by, 

4> 2 SS - {2A VV (1 - a) - 1 + ap)4> ss - ap = (29) 



while the normalized mean fitness for the population 
replicating with the random sexual replication strategy 
is given by, 

2A vv {l-a) 2 <p 2 rs -(<j> rs + l)((f> rs -ap)(<p rs (l-2a)+ap) = 

(30) 

The central result of this subsection is that 4> ss > (f> rs > 
A)$ a , except when p = 0, p = 1, a = 0, or a = 1, in which 

case 4> ss = (j) rs = 4>a- 
(26) -\Y e w iH prove this result in two steps: First we will 
prove that <j) ss = 4> rs = <f> a for p = 0, p = 1, a = 0, or 
a = 1. We will then prove that (p ss > <p rs > <p a as long 
as a,p € (0, 1). 

1. Proof that (j> ss = 4> rs = 4> a at p = 0,1 and/or a = 0, 1 

When p = 0, <j) ss , <p rs are obtained by solving, 

(p ss {4> ss + 1 - (1 - a)s) = 
# a (*(l - a) 2 - (1 - 2a){<j> rs + 1)) = (31) 

which have the solutions 4> ss = 4> sr = 0. We choose these 
solutions, because they are the ones that are physical. 
When p = 1, (f> ss , <p rs are obtained by solving, 

(<j> ss - 1)(4> SS + a) = 
{<j> sr ~l){<f sr {l-2a) + a 2 ) -0 (32) 

so that (j> ss = 4> sr = 1. 
When a = 0, we obtain, 

4> ss ((t>ss - (2A VV - 1)) = 
<j> 2 rs (cp rs - (2A VV - 1)) = (33) 

These two equations both admit the solutions and 
2A VV — 1. Since <f) ss = 4> rs = 1 for p = 1, and 
(p ss = 4> rs = for p = 0, by continuity it follows that 
4> ss = 4> rs = 2A VV - 1 for p e [p cri t,l], and for 

p e [o,p cri t]. 

When a = 1, we obtain, 

(<i>„ - p)(<I>,8 + 1) = 

(4> rs - P ) 2 (^ rs + l) = (34) 

so that 4> ss — <ft rs = p. 

Note that in all cases, we have (f> ss = <f) rs = <j> a , as we 
wished to show. 



2. Proof that <f) 3S > <f) rs > <f> a when a,p € (0, 1) 

When a,p E (0, 1), we claim that Eq. (30) has a so- 
lution in (ap,p). By continuity, we expect that this so- 
lution is the value of cf) rs as a function of p, since it is 
the solution that is consistent with <j> rs = at p = and 
(f) rs = 1 at p = 1. 

When 4> rs = ap, Eq. (30) evaluates to 2A vv (l — 
a) 2 a 2 p 2 > 0. When 4> rs — p, Eq. (30) evaluates to 
-(1 - a) V(l - p)(l - *(1 - P)) < 0. 
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By the Intermediate Value Theorem, it follows that 
Eq. (30) has a solution in the interval (ap,p). This of 
course shows that <fi rs > ap. 

We now claim that (f) rs ^ 2A VV — 1. For, if (fi T s — 
2A VV — 1, then from Eq. (30) we have, 



{(j>rs-p) 2 = 



(35) 



which implies that <j) rs — p. But this means that 2p + 
s(l-p) 2 -l =p=> (l-p)(l-s(l-p)) = =>p= 1 
since p € (0, 1) by assumption. Therefore, <f) rs ^ 2A VV — 
1, as claimed. 

If we can now show that there exists a p € (0, 1) such 
that 4> rs > 2A VV — 1, then we will have proven that 
phi rs > 2A VV — 1 for all p e (0,1). For otherwise, by 
the Intermediate Value Theorem we would be able to 
find ape (0, 1) such that (f) rs = 2A VV — 1 

Now, p cr it is defined by the equation 2p crit + s(l — 
Pent) 2 - 1 = apent, so that p crit = 1 a = 1 and 
Peru = => s = 1. Since we are assuming a S (0,1), 
we either have s < 1 or s = 1. If s < 1, then p crit g 
(0, 1), so since 2A VV — 1 = ap at p = p cr it, we have that 
rs > 2A VV - 1 at p = p cr it e (0,1), thereby proving 
that cf) rs > 2A VV - 1 for p e (0,1)- Taking the limit 
s — > 1 gives that rs > 2A TO — 1 for p e (0, 1), which of 
course implies that (firs > 2A VV — 1 for p £ (0, 1), since 
rs # 2A TO - 1 forpe (0,1). 

We have now shown that cfi rs > ap,2A vv — 1 when 
ci,pe (0, 1), and so <fi rs > <fia for a,p G (0, 1). We now 
turn to proving that <fi ss > <fi rs for a,p € (0, 1). 

If we define f2 = 2A VV (1 — a) — 1 + ap, then </>g S = 
Q(fiss + cep. If (fi ss = (fi S r for some p G (0, 1), then defining 
= S;S = </> sr , we obtain from Eq. (30) that, 

2A VV (1 - a) 2 (n<fi + ap) 

= <fi(n + 1 - ap){(fi{\ - 2a) + ap) 

= + 1 -ap)[{{\ -2a)n + ap)(fi+ (1 -2a)ap] 

=> [2A vv (l - a) 2 fl -(Q + l- ap){{\ - 2a)Q + ap)]<fi 

= (Q + 1 - ap)(l - 2a)ap - 2A VV (1 - afap (36) 

Now, ft + 1 — ap = 2A VV (1 — a), so that, 

[n-p](fi=-ap (37) 



So, multiplying Eq. (29) by (ft — p) 2 gives, 

+ - p)ap - (ft - pfap = 
(l-a)(l-p)(l-s(l-p)) 



2 2 

a p 



= 




(38) 



which is impossible for a,p e (0, 1). Therefore, 
for p e (0, 1). 

Now, let us look at d<fi ss /dp and d<fi rs /dp at p = 1. For 
0.,., we have, 



while for 



dp 
we have, 



1 + a 



^)- = 2 



(39) 



(40) 



Therefore, for p near 1, we have, 

2 



= 1 - 



1 + a 
b rs = 1-2(1- p) 



(1-P) 



(41) 



so since 2 > 2/(1 + a), it follows that <fi ss > cfi rs forp < 1, 
but p close to 1. But then (f> ss > <fi rs for all p e (0, 1), for 
otherwise there exists a p G (0, 1) for which <fi ss < <fi rs , 
which implies by the Intermediate Value Theorem that 
there exists a p € (0, 1) such that <fi ss = <fi rs => < ^=. 

We have now established that <fi ss > cfi sr > (fi a for 
a,p e (0, 1), as we set out to prove. 



B. Case 2: n vv /(-yp) > 

The central result of this subsection is that <fi ss > <fi rs 
for a,p G (0, 1), with equality only occurring when p = 0, 
p = 1, a = 0, or a = 1. 

Suppose we can show that <fi ss ^ 4> rs for a,p € (0,1), 
independently of the value of K vv /(-fp). Then if (fi ss < <fi rs 
for some a,p G (0, 1) and for some value of K vv /{^p) > 
0, it follows from continuity and the Intermediate Value 
Theorem that <fi ss — (firs for some value of n vv /(jp) > 0, 
since <fi ss > (firs for a,p£ (0, 1) when n vv /(-fp) = 0. This 
is a contradiction, and so we must have that cfi ss > <fi rs 
for all a,p E (0, 1) and n vv /(jp) > 0. 

Therefore, to prove the central result of this subsection, 
we will show that <fi ss — <fi rs for p = 0, 1 and/or a = 0, 1, 
and that cfi ss ^ <fi rs f° r a,P <= (0, 1). 

Working again with the definition A — K vv /(2jp), we 
have from Eq. (14) that, 



((fi 2 s - (2A vv (l - a) - 1 + ap)(fi ss - apf 
(0ss + l) 2 ^ss + a) 2 



= Hss (42) 



As A increases from to oo, we expect (fi ss to decrease 
from its maximal value down to 0. In this regime, we 
would like to determine the sign of the expression 4> 2 S — 
{2A VV {1 - a) - 1 + ap)<fi ss - ap. 

Consider the polynomial x 2 — Ax — B, which has the 
roots (l/2)[A±VA 2 + 4:B}. If B > 0, then the "+" root, 
denoted r + , is positive, and the "-" root, denoted r_, is 
negative. Therefore, writing x 2 — Ax — B = (x — r_)(x — 
r + ), we may note that when x e [0, r + ], x — r_ > and 
x — r + < 0, so that x 2 — Ax — B < 0. 

The positive root of <fi 2 s — (2A vv (l — a) — l + ap) <fi ss — ap 
is the value of (fi ss when A — 0. Since (fi ss then decreases 
to as A increases, it follows from our analysis that <fi 2 s — 
(2A vv (l — a) — 1 + ap)<fi ss — ap is zero or negative for 
A > 0. Therefore, taking the square roots of both sides 
of Eq. (42) gives, 



ap+ (2A vv (l - a) - 1 + ap)(fi s 
(<j) ss + l)((j)ss + a) 



= ^M7s (43) 



If, 



for 



lui a given a,p G [0, 1] and A > 0, then we 
may insert the expression above into Eq. (27) defining 
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4> rs . If we set <j> — <fiss — <firs, then we obtain, after some 
manipulation, 

= [0(1 - 2a) -a]x 

[<P 2 - (2A vv (l - a) - 1 + ap)<P - ap] 2 

-2[4> + l] x 

[(1 - 2a)ct) 2 -a 2 (l- p)<f> + a 2 p] x 

[(f) 2 — (2A VV (1 — a) — 1 + ap)4> — ap] 

+2A vv (l - a) 2 4> 2 [4> 2 - {2A VV {1 - a) - 1 + ap)4> - ap] 

-2A vv {l-a) 2 <p 2 [<P+l][c^ + a] 

+ $+l] 2 $ + a][(f>-ap][(f>(l-2a)+ap] = (44) 

This expression may be simplified (with the aid of a 
symbolic math package if necessary) to give, 

= a(l-a) 2 (l-p)^(l_0 s ) x 

[r 2 p 3 + 3r(l - r)p 2 + (2(1 - r) 2 + r 2 )p + r(l - r)] 

(45) 

which implies that either a = 0, 1, p = 1, or <f> = 0, 1. 
Since </> = is equivalent to p = when a > 0, and since 
= 1 is equivalent to p = 1 and A = 0, we obtain that 
(f) ss = <p rs only when a = 0, 1 or p = 0, 1. If a,p € (0, 1), 
then cj)gg ^ <f) rs . 

Therefore, we have proven that <j> ss — 4> rs when p = 0, 
p = 1, a = 0, or a = 1, and SS > </) rs for a,p € (0, 1), 
independently of the value of n vv /(-fp). 

V. DISCUSSION 

The two key results of this paper are that a sexual 
population employing a random mating strategy will out- 
compete an asexual population when the cost for sex is 
negligible, and that a sexual population using a selective 
mating strategy will outcompete a sexual population us- 
ing a random mating strategy. The only exceptions are 
the boundary cases p = 0, 1 and a = 0, 1. However, even 
here, the mean fitnesses of the two sexual strategies are 
identical. Furthermore, when the cost for sex is negligi- 
ble, then the mean fitnesses of the sexual strategies are 
identical to that of the asexual strategy. 

That random mating provides a selective advantage 
over asexual replication is an interesting result, because 
the strategy can lead to the formation of diploids with 
completely defective genomes. Presumably, however, the 
fitness benefit provided by the formation of diploids with 
two functional chromosomes outweighs the fitness cost 
associated with the formation of diploids with two de- 
fective chromosomes, leading to an overall advantage for 
the strategy. Nevertheless, because the selective mating 
strategy does not produce genomes with defective chro- 
mosomes, this strategy has an advantage over the random 
mating strategy. 

When the cost for sex is negligible, an analysis of the 
mean fitnesses <f> a , (f> ss , and <j) rs near p = 1 yields some 
interesting results. We have <t> a — 2A VV — 1 = 1— 2(1 — 



p) + s(l — p) 2 . We have also shown that when p is close to 
1, then to first-order in l-pwe have <j) rs = 1 — 2(1 — p) 
and <j) ss = 1 — (2/(1 + a))(l — p). Therefore, we may 
note that the random mating strategy and the asexual 
strategy are identical to first-order in 1 — p, while the 
selective mating strategy already outcompetes both the 
random and asexual strategies. 

However, from Eq. (30) we obtain, 

(^W=2<« + 2 (r ^) (46) 
and so, to second-order in 1 — p we have that, 

<j> r „ = 1 - 2(1 -p) + [s + 2(-^-) 2 ](l - p) 2 (47) 

1 — a 

Comparing the second-order expression for <j> rs to <f> a , we 
see that <fi rs exceeds <p a by 2(a/(l — a)) 2 (l — p) 2 when p 
is close to 1. 

Figure 3 shows a plot of the two sexual replication 
strategies when there is no cost for sex and the asexual 
strategy. 

When there is a cost for sex, then when p = 1 the 
asexual population outcompetes both sexual populations. 
This makes sense, for when p — 1 replication is error-free, 
and hence at steady-state the asexual population consists 
only of the wild-type. In this case, genetic recombination 
will not improve fitness, since there are no defective chro- 
mosomes in the population to begin with. 

However, because sexual replication when there is no 
cost for sex will outcompete an asexual replication when 
p € (0, 1), it follows that if the cost for sex is sufficiently 
low, then below a certain value of p a sexual population 
will outcompete the asexual population. Presumably, the 
higher the cost for sex, the smaller p must be before the 
selective advantage for sexual replication is sufficiently 
large to outweigh the cost. 

This behavior only persists up to a maximal cost for 
sex, beyond which asexual replication outcompetes sex- 
ual replication at all replication fidelities. The reason for 
this is that the sexual and asexual mean fitnesses con- 
verge to as p — > 0. As a result, once the cost for sex 
is sufficiently high, the fitness advantage of the sexual 
strategy for the values of p where the sexual strategy can 
outcompete the asexual strategy is too small to overcome 
the cost for sex. 

Indeed, from Eq. (14) it may be shown that 
(d(j) ss /dp) p= o = when K vv /(-yp) > 0. Since 
(d<j) a /dp)p=o = a, it follows for a e (0, 1) that <j) ss < 4> a 
in a neighborhood of p sufficiently close to 0. Since 
4>rs _ 4>ss, the same condition holds for c/) rs as well. How- 
ever, if any of the sexual populations can outcompete the 
asexual population when p is sufficiently large, then there 
must exist another value of p below which asexual repli- 
cation outcompetes sexual replication, and above which 
sexual replication outcompetes asexual replication. 

The complete picture is then one where, for a non-zero 
cost for sex, the asexual population outcompetes a sexual 
population above a certain replication fidelity. Below this 
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FIG. 1: Regimes where the asexual, selective sexual, and 
random sexual replication strategies are respectively advanta- 
geous, as a function of p and K, vv /("/p). The region below the 
dotted curve is the region where both the random and selec- 
tive mating strategies outcompete the asexual strategy. The 
region below the solid curve and above the dotted curve is the 
region where only the selective mating strategy outcompetes 
the asexual strategy. The region above the solid curve is the 
region where the asexual strategy dominates. Note that as 
the cost for sex increases, the regions where the sexual strate- 
gies are advantageous shrink and eventually disappear. The 
random mating strategy ceases to be advantageous at a lower 
cost than the selective strategy. The parameters chosen are 
a = r= 1/2. 

replication fidelity, the sexual population outcompetes 
the asexual population. Finally, once the replication fi- 
delity becomes sufficiently low, the asexual population 
again outcompetes the sexual population. 

Based on this analysis, we expect that, as the cost 
for sex increases from to oo, the region of replication 
fidelities where sexual replication outcompetes asexual 
replication starts at (0,1), gradually shrinks, and even- 
tually disappears once the cost for sex crosses a thresh- 
old value. Of course, because the random mating strat- 
egy has a lower fitness than the selective mating strategy, 
the region of replication fidelities where the random mat- 
ing strategy outcompetes asexual replication is a subset 
of the region of replication fidelities where the selective 
mating strategy outcompetes asexual replication. 

Figure 4 shows a plot of the asexual and sexual mean 
fitnesses when there is a non-zero cost for sex. Figure 
5 shows a plot of the regions for the selective advan- 
tages of the various replication strategies as a function of 

Kvv/{lP)- 

Although we have shown that a selective mating strat- 
egy will outcompete a random strategy, this analysis is 
based on the assumption that a selective mating strategy 
is not inherently slower than a random mating strategy. 
As discussed in the Introduction, this is in general not 
true, since there is a time cost associated with determin- 
ing the genome of a potential haploid mate. This time 
cost leads to an additional fitness cost associated with a 
selective mating strategy that is not incurred by a ran- 



dom mating strategy. When the replication fidelity is 
either sufficiently high or sufficiently low, then the ben- 
efits of a selective mating strategy are not sufficient to 
overcome the fitness costs, and so the random mating 
strategy dominates. However, at intermediate replica- 
tion fidelities, the selective mating strategy may indeed 
outcompete the random mating strategy, assuming that 
the fitness cost associated with selective mating is not too 
high. Once the fitness cost associated with selective mat- 
ing becomes sufficiently high, then the random mating 
strategy may outcompete the selective mating strategy 
at all replication fidelities. 

In summary, when the time cost for a selective mating 
strategy is taken into account, an analysis of the regimes 
where the selective and random mating strategies are ex- 
pected to be respectively dominant may produce a curve 
that is analogous to the ones in Figure 5. 



VI. CONCLUSIONS AND FUTURE RESEARCH 

This paper developed a set of simplified models de- 
scribing asexual and sexual replication in a unicellu- 
lar population consisting of two-chromosomed, diploid 
genomes. We considered two types of sexual replication 
strategies: A selective mating strategy, where only viable 
haploids are allowed to fuse, and a random mating strat- 
egy, where all haploids are allowed to participate in the 
replication process. We assumed that haploid fusion was 
a second-order rate process. 

We found that, when the cost for sex is negligible, 
both the selective and random mating strategies lead to 
a greater mean fitness than asexual replication. Never- 
theless, we found that the selective mating strategy has 
a higher mean fitness than the random mating strategy, 
as long as the additional fitness penalty associated with 
a selective mating strategy is negligible. 

Further analysis suggested that sexual replication for 
both the selective and random mating strategies is fa- 
vored at intermediate mutation rates, and when the cost 
for sex is sufficiently low. Once the cost for sex becomes 
sufficiently high, the selective advantage for sexual repli- 
cation disappears entirely. 

The results of this paper therefore suggest that sex is 
favored in slowly replicating organisms and high popu- 
lation densities. While this is consistent with previous 
studies [9l [lOj HH E2], what is interesting is that this 
result holds even with a random mating strategy. 

For future work, we will consider models where it is 
not necessarily true that n uu = 0, but rather we will al- 
low for positive values of k uu . This will make our model 
more consistent with quasispecies models that often as- 
sume a small, but positive growth rate for organisms with 
defective genomes. 

Furthermore, although our model provides results that 
are broadly consistent with actual organismal behavior, 
we have nevertheless worked with a highly simplified sex- 
ual replication model that does not exactly correspond to 
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the sexual replication pathway employed by unicellular 
organisms. 

In our model, a mature diploid divides into two hap- 
loids, the haploids fuse, and then the resulting diploid di- 
vides into two cells. Thus, the diploid mitosis occurs after 
haploid fusion. By contrast, in Saccharomyces cerevisiae 
(Baker's yeast), for example, a mature diploid divides 
into two diploids, which then divide into four haploids. 
The haploids then fuse, and the resulting diploids grow 
to maturity to begin the process again. Thus, for yeast, 
diploid mitosis occurs before haploid fusion. We suspect 
that the relative positions of the diploid mitosis and hap- 
loid fusion stages in a sexual replication cycle can affect 
the mutation-selection balance, and so we would like to 
develop more realistic models corresponding to the sexual 
replication pathways in actual unicellular organisms. We 
should emphasize, however, that we believe that the cen- 
tral results of this paper will hold even when we consider 
a somewhat different sexual replication pathway. 



Finally, one of the simplifying assumptions made in 
this paper is that the organisms have genomes consist- 
ing of two chromosomes. For future research, we plan to 
study genomes consisting of arbitrary numbers of chro- 
mosomes. Within this context, we plan to consider more 
complex fitness landscapes and the role of intra-genomic 
recombination. 
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